Signatures of Discontinuity in the Exchange-Correlation Energy Functional Derived 
from the Subband Electronic Structure of Semiconductor Quantum Wells 
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The discontinuous character of the exact exchange-correlation (xc) energy functional of Density 
Functional Theory is shown to arise naturally in the subband spectra of semiconductor quantum 
wells. Using an ab-initio xc functional, including exchange exactly and correlation in an exact par- 
tial way, a discontinuity appears in the xc potential, each time a subband becomes slightly occupied. 
Exchange and correlation give opposite contributions to the discontinuity, with correlation overcom- 
ing exchange. The jump in the intersubband energy is in excellent agreement with experimental 
data. 
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Density Functional Theory (DFT) has become the 
standard calculational tool in physics and quantum 
chemistry for the study of atomic, molecular, and solid 
state systems. The theory is based in the Hohenberg- 
Kohn theorems that place the ground-state electron 
density as the basic variable and provides a variational 
principle for its calculation. Kohn and Sham (KS) 
showed how the problem of variational minimization for 
the density could be exactly mapped to one of non- 
interacting particles in an effective potential, which con- 
tains only one non-trivial component: the exchange- 
correlation (xc) contribution[IJ. DFT, however, gives 
no clue on how to proceed for its practical calculation. 
Naturally a lot of attention has been devoted to the 
development of better xc energy functionals; KS intro- 
duced the highly successful Local Density Approximation 
(LDA), which is widely employed nowadays, along with 
the improvements born from it (such as the Generalized 
Gradient Approximation or GGA, meta-GGA, etc.)[2]. 
The work described here is motivated by this funda- 
mental need of better approximations for the xc energy 
functional 0|, using as "laboratory" to test the accuracy 
of the approximations the subband electronic structure of 
the quasi two-dimensional electron gases (2DEG) formed 
at the interface between two dissimilar semiconductors, 
such as GaAs and AlGaAs. In this Letter we show that 
at the one-subband — > two-subband quantum well (QW) 
transition (15 — > 25), the xc potential behaves discon- 
tinuously, with exchange and correlation giving opposite 
contributions (i.e., competing) to the discontinuity. The 
intersubband energy, which also jumps abruptly at the 
transition, is in excellent agreement with experiments. 

Our model system is a semiconductor modulation- 
doped QW grown epitaxially, as shown in the upper in- 
set of Fig. 1. Assuming translational symmetry in the 
(x — y) plane (area ^4), and proposing accordingly a solu- 
tion of the type <^ik<r(p, z) — e lk 'P£![(z)/\/A, with p the 
in-plane coordinate, the zero temperature ground-state 
electron density can be obtained by solving a set of effec- 
tive one-dimensional KS equations of the form: 



where effective atomic units have been used. £f (z) is the 
wavefunction for electrons in subband i (i — 1, 2, ...), 
spin a (a =t, I), and eigenvalue e\. The local (mul- 
tiplicative) KS potential Vj£ s (z) is the sum of several 
terms: Vg 8 {z) = V ext (z) + V H {z) + V° c (z). V ext {z) is 
given by the epitaxial potential plus an external elec- 
tric field. Vh{z) is the Hartree potential. Within DFT, 
V xc (z) = A~ 1 5E XC / 5n a (z). Departing from the main 
stream in most applications of KS-DFT, our xc energy 
functional is an explicit functional of the whole set of 
ef's and £f's, E xc = E xc [{ef , £f}], but an implicit (in 
general unknown) functional of the spin-resolved den- 
sity n a (z). The zero-temperature 3D electron density 
is n a (z) = Yl e 7<n(f x — £ f )|£f ( z )| 2 /2tt, with (j, the chemi- 
cal potential. By assuming a paramagnetic situation, we 
drop the spin index a from this point. 

The xc energy functional E xc for our 2DEG has been 
generated by Gorling-Levy (GL) perturbation theory 
where the correlation energy is expanded in a series [3], 



E xc [{e i) Z i }] = E x [{£ i ,Z i }} 
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which is truncated at its leading contribution n = 2. E x 
in equation above is the exact exchange energy, known 
explicitly as a Fock integral of the KS occupied orbitals. 
Its multi-subband explicit expression is given by Eq.(38) 

in Ref.Q (denoted as I below). The terms E^ L ^ can 
be found explicitly Q. At difference with E x , however, 
they depend both on occupied and an infinite number 
of unoccupied subbands. Their numerical evaluation is 
in consequence rather expensive. The explicit expression 

GL (2) 

for E" y ' for our semiconductor QW system is given 
by Eqs.(32) and (33) of I. A similar correlation energy 
functional has been used in Refs. and for the case 
of atoms, with mixed results. As shown here, Eq.(2) 
seems much more promising for the ^Z)SG0|. 

The next non-trivial problem is the evaluation of 
V xc (z), given the already quoted implicit dependence of 
E xc on n(z) in Eq. (2). The procedure for dealing with 
implicit functionals relies on the use of the chain rule for 
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Figure 1: Ground-state —* first-excited intersubband energy 
£12 as a function of F. Negative (positive) values of T corre- 
spond to the IS (2S) regime. Upper inset: schematic view of 
our model for the modulation-doped QW. A charge-transfer 
field along z is induced by a distant metallic plate (thin verti- 
cal line at right). If the metallic plate is positively (negatively) 
charged, more (less) electrons are transferred towards the well 
from an ionized donor impurities region (thick vertical stripe 
at left), which acts as a particle reservoir fixing the chemical 
potential fi. Lower inset: V xc (z) in the same three different 
approximations as for £12 (r = 0.5 meV) . The QW extends 
from z = 50 A to z = 295 A. 



functional derivatives as follows 



V xc {z) 



1 5E XC 
A Sn(z) 



SE XC 8V K s{z') 
SV K s(z') 6n(z) 



dz'. (3) 



To proceed with the calculation of V xc (z) directly from 
Eq.(3), we use a numerical method devised and explained 
in detail in I. Eqs.(l)-(3) should be iterated until full self- 
consistency is achieved. 

We present in Fig. 1 the intersubband energy spac- 
ing £12 = £2 — £1, as a function of Y = ^ — £ 2 , for 
three approximations for E xc : LDA, exact-exchange plus 
LDA for correlation [X(X)-C(LDA)], and X(X) plus par- 
tial exact-correlation (E c '), denoted as X(X)-C(pX). 
Note that while the resulting £12 are quite similar in the 
three approximations in the whole 15 regime, and in the 
25 regime with sizable occupation of the second sub- 
band, noticeable differences arise in the limit of small 
second subband occupancy. Starting with the X(X)- 
C(LDA) approximation, £12 shows an exchange-driven 
abrupt positive jump at the 15 — > 250. The inclusion 
of C(pX) overcomes the X(X)-C(LDA) positive jump, re- 
sulting now in a negative jump in £12, until it levels with 
the other results at finite second subband fillings. LDA 



is in between the two kind of discontinuities at T = 0, 
showing only a discontinuity in the derivative. The lower 
inset shows the corresponding V xc (z) in the relevant QW 
region. V xc (z) as resulting from X(X)-C(LDA) builds a 
barrier just where most of the weight of £2(2) is concen- 
trated, pushing £2 upwards under small occupancy of the 
second subband, and leaving S\ more or less unaltered; 
this explains the discontinuous positive behavior of £12. 
It is also physically reasonable: by blocking the occupa- 
tion of the second subband, the exchange energy is opti- 
mized, as the intrasubband exchange is larger than the 
intersubband exchange. The V xc (z) in the X(X)-C(pX) 
behaves just in the opposite way: it develops a deep well 
just after the transition, inducing an abrupt decrease of 
£2 at more or less constant £1. This explains the abrupt 
decrease of £12 in this case. The behavior has again a sim- 
ple physical explanation: by inducing the occupancy of 
the second subband, V xc (z) promotes a spatial separation 
between electrons in both subbands, decreasing correla- 
tion and its associated repulsive energy. With respect 
to the V xc (z) resulting from LDA, it shows the expected 
smooth and continuous behavior at the transition. 

Besides these fully numerical results, we provide below 
an analytical derivation of the results of Fig. 1, for T ~ 0. 
In the limit |T| — > 0, E xc can be approximated as 



e£>=p£> + tq 



(«) 

xc ' 



with a = 15,25. Here, P^ C S) = E XC (T 

p( 2S ) _ TP (y 

1 xc — J -/ Xc y± 



0+), Q X \ S) = dE xc /dT\ -, Q [ x 7> = 
dE xc /dT\ +. By inspection of the explicit expressions for 

GL(2) (IS) 

E x and E c given in I, it is concluded that P x 
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P xc , and that Qxc 
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In words, for a 



fixed set of ej^'s and £j(z)'s, the xc functional is contin- 
uous at the 15 — > 25, but its derivative is discontinuous. 
The explicit expressions for P xc , Q^ 1 ^ and Q X 2 C S ^ are not 
needed for the present derivation. Inserting Eq.(4) in 
Eq.(3) we obtain, 



SP XC 



5V KS {z') 



X - 1 (z' 1 z)dz'. (5) 



Here, we have defined XaH z i z ') = &Vks(z) / '6n a (z'), 
and used the result ST j 'SVks(z) = —8s2/8Vks(z) = 
— \&(z)\ 2 , obtained by first-order perturbation theory. 
We have also neglected a term lineal in T, which becomes 
arbitrarily small in the limit T — ► 0. In the limit r — > + , 
the density response function^ becomes 



X2s(z, z') = xis(z, z') - Mz)W)\ fa- 



(6) 



Its inverse can be calculated by using the Sherman- 
Morrison technique|10|. We obtain 

X2S 1 (2, A = Xri(*. z ') + D(z)D(z')/[n(l + A)], (7) 

where D(z) — J Xis( z i x )\&(x)\ 2 dx, and A = 
— 7T _1 J D(a;)|^2(^)| 2 ^- As we can see from Eq.(7), 
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X2s( z ' z ') m discontinuous at the 15 — > 25 transition, 
such as it is X2s{z,z') of Eq.(6). Using now Eq.(5) for 
a = 15, 2S and exploiting the explicit expression for 
X2S ( z > z ') as given by Eq.(7), we arrive at the result 

AV xc (z) = v£ s \z)-V£ s \z) = C xc D(z)/(l + \), (8) 

with C xc = (v£ S) ) 2 /ir-AQ xc , (0} l = J dx\^(x)\ 2 0(x) 

and AQ XC = Q x 2 f ^ — Q^fK Eq.(8) is an important re- 
sult, which shows explicitly how the functional depen- 
dence on z of the xc potential changes discontinuously 
at the 15 — * 25 transition. For an arbitrary subband 
transition NS — * (N + 1)5 it can be shown that the re- 
sult of Eq.(8) is still valid, with the replacements 1 — > N 
and 2 — > N + 1 in Eqs.(4)-(8). Equation 8 follows rig- 
orously from Eq.(4), and then it is important to discuss 
its validity. In writing Eq,(4) we have implicitly made a 
"frozen" assumption, by taking the same set of wavefunc- 
tions and energies e% both for r — > 0~ and T — > + . 
In the language of the numerical self-consistent calcula- 
tions leading to Fig. 1, it is as if the results for the 15 
case in the limit r — > 0~ were extrapolated to the 25 
r — > + case, by performing a single iteration loop to- 
wards full self-consistency. This "frozen" result for the xc 
energy functional is illustrated in Fig. 2 by the straight 
lines, separating exchange (P x + TQ") from correlation 
(Pc + TQ"). Besides, if the approximation for E xc is such 
that AV xc (z) 7^ (that is, if C xc ^ 0), it is clear that 
the self-consistent iteration loop will lead to a further 
discontinuity in the xc energy functional itself right at 
the 15 — * 25 transition (r = 0), once convergence has 
been reached. These fully self-consistent results corre- 

GL ( 2") 

spond to the thick full (E x ) and dashed (E c ') lines 
in Fig. 2. On the other side, if the approximation for 
E xc is such that C xc = 0, no discontinuity of the type 
of Eq.(8) exists for V xc (z) before or after self-consistency 
is achieved, which in turn implies the continuity of E xc . 
Our E xc = E xc [{ef, £f}] is such that C xc ^ 0; local 
(LDA) approximations for E xc yield C xc = 0. 

For further insight, the constant C xc can be separated 

conveniently in its exchange and correlation components 

(is) 

C n = (Vn — AQ n , with n = x,c. The analysis of 

the results shown in Fig. 2 leads to the conclusion that 
C x < and C c > 0, considering that D(z) < (see inset 
in Fig. 2). Qualitatively, this happens in the following 

way: a) in the exchange case, (Vx )a is large and neg- 
ative (see inset Fig. 2), while —AQ X is a relatively small 
positive magnitude, resulting in C x < 0; b) in the corre- 

lation case, (V c /2 is a very small quantity (see inset 
Fig. 2), while — AQ C is a relatively large positive number, 
yielding a C c > 0. The net result is that correlation over- 
comes exchange (C xc = C x + C c > 0), and AV xc {z) has 
a negative contribution resulting in the right attractive 
well shown in the lower inset of Fig. 1. It is important to 
note that this overcoming of correlation on exchange hap- 
pens even when \E C y '\ <§; \E X \ (Fig. 2). However, as 
Eq.(8) clearly shows, not only the magnitude of the xc 



-0.5 




r = u-e 2 [meV] 



Figure 2: Full lines: E x in the "frozen" (thin line) and self- 
consistent (thick line) approaches; dashed lines: E^ L ^ in the 
"frozen" (thin line) and self-consistent (thick line) approaches. 
Note the cut in the vertical scale. Inset: exchange and cor- 
relation potentials (full lines) and function D(z) (dotted line, 
arbitrary units) for V = — 0.01 meV. 
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energy functional matters (represented by the (V xc )2 
contribution), but also the respective derivatives (repre- 
sented by the AQ XC term). This exemplifies quite vividly 
the potential danger of neglecting correlation against ex- 
change at the threshold of a subband transition, under 
the argument that the correlation energy is much smaller 
than the exchange energy. 

Is there any experimental evidence of this type of dis- 
continuities in semiconductor QW's? The answer is yes. 
We reproduce as an inset in the upper panel of Fig. 3 the 
experimental values for the subband densities n\ and 112, 
plotted as a function of the Fermi level measured from 
the top of the valence band[llj. The data have been ob- 
tained from a quantitative analysis of photoluminescence 
line shapes. It is seen that when the Fermi level touches 
the bottom of the second subband, the electron density 
n<2 jumps from zero to a finite value ranging from 3x 
10 10 to 8 x 10 10 cm~ 2 , depending on temperature[ll|. 
The electron density of the lowest subband, in contrast, 
increases slightly but smoothly with voltage, suggesting 
that the external electric field couples essentially to the 
15 occupation, which seems quite plausible as it is the 
subband with the largest occupation. This motivates us 
to redraw the results for £12 of Fig. 1 as a function of 
PL— e\, in the lower panel of Fig. 3. Clearly the theoretical 
n\, 122, and £12 vs. ;U — E\ agree quite well with the exper- 
imental data, both qualitatively and quantitatively. For 
instance, the theoretical value for the decreasing jump 
in £12 at fx — £1 ~ 25 meV, amounts to about 3.3 meV, 
in excellent agreement with the 3.5 meV jump estimated 
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Figure 3: Thick full lines: Ground and first-excited subband 
densities n\ and 712 (upper panel) and intersubband transition 
energy £12 (lower panel) as a function of /U — £1. Dashed lines: 
n2 and £12 corresponding to a continuous filling of the first- 
excited subband, as in Fig. 1. Insets: experimental data from 
Ref.fllll. 



from experiment. 

The results presented in this work are intimately re- 
lated to the issue of the derivative discontinuity of ensem- 
ble DFT(i^|. Among the many important consequences 
derived from this extension of DFT to fractional particle 
number, maybe the most important one in the field of 
solid state physics has been the realization that the true 
gap of semiconductors is not given by the KS one-particle 
gap[3|. Instead, the true gap is given by the sum of the 



KS gap, plus the so-called xc discontinuity, A xc . Contin- 
uum approximations to xc energy functionals, including 
all currently widely used LDA's and GGA's, fail to pro- 
duce the correct value for A xc , resulting in an important 
underestimation of the fundamental band-gap of most 
semiconductors and insulators. In a very recent work, 
Griming et al have clarified the theoretical situation, ob- 
taining very good agreement between experimental and 
theoretical band gaps of Si, LiF, and ArflJ, by using an 
orbital based correlation energy functional correspond- 
ing to a dynamical screening of the Coulomb interaction 
(G W approximation) . It was found that A xc contributes 
as much as 30% — 50% to the energy gap. Their results, 
for a different type of systems, are fully consistent with 
ours. 

In conclusion, the intrinsic discontinuity of the xc 
energy functional has been obtained entirely within a 
DFT framework, for a realistic system. We have shown 
that the xc energy functional generated by second-order 
Gorling-Levy perturbation theory for the 2DEG has 
many of the properties of the exact functional. The main 
finding is that at the IS — > 2S transition, the KS xc po- 
tential and the associated intersubband transition ener- 
gies behave discontinuously, with x and c giving oppo- 
site contributions (i.e. competing) to the discontinuity, 
and with correlation overcoming exchange. Very good 
qualitative and quantitative agreement is obtained with 
experiments. 
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